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Abstract 

High-value transactions between banks in Australia are settled in the Reserve Bank Information and Trans- 
fer System (RITS) administered by the Reserve Bank of Australia. RITS operates on a real-time gross 
settlement (RTGS) basis and settles payments and transfers sourced from the SWIFT payment delivery 
system, the Austraclear securities settlement system, and the interbank transactions entered directly into 
RITS. In this paper, we analyse a dataset received from the Reserve Bank of Australia that includes all 
interbank transactions settled in RITS on an RTGS basis during five consecutive weekdays from 19 February 
2007 inclusive, a week of relatively quiescent market conditions. The source, destination, and value of each 
transaction are known, which allows us to separate overnight loans from other transactions (nonloans) and 
reconstruct monetary flows between banks for every day in our sample. We conduct a novel analysis of 
the flow stability and examine the connection between loan and nonloan flows. Our aim is to understand 
the underlying causal mechanism connecting loan and nonloan flows. We find that the imbalances in the 
banks' exchange settlement funds resulting from the daily flows of nonloan transactions are almost exactly 
counterbalanced by the flows of overnight loans. The correlation coefficient between loan and nonloan im- 
balances is about —0.9 on most days. Some flows that persist over two consecutive days can be highly 
variable, but overall the flows are moderately stable in value. The nonloan network is characterised by a 
large fraction of persistent flows, whereas only half of the flows persist over any two consecutive days in the 
loan network. Moreover, we observe an unusual degree of coherence between persistent loan flow values on 
Tuesday and Wednesday. We probe static topological properties of the Australian interbank network and 
flnd them consistent with those observed in other countries. 
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1. Introduction 

Financial systems are characterised by a complex and dynamic network of relationships between multiple 
agents. Network analysis offers a powerful way to describe and understand the structure and evolution of 
these relationships; background information can be found in and 0]. The network structure plays an 

important role in determining system stability in response to the spread of contagion, such as epidemics in 
populations or liquidity stress in financial systems. The importance of network studies in assessing stability 
and systemic risk has been emphasised in Q in the context of integrating economic theory and complex 
systems research. Liquidity stress is of special interest in banking networks. The topology of a banking 
network is recognised as one of the key factors in system stability against external shocks and systemic 
risks 0. In this respect, financial networks resemble ecological networks. Ecological networks demonstrate 
robustness against shocks by virtue of their continued survival and their network properties are thought to 
make them more resilient against disturbances Often they are disassortative in the sense that highly 
connected nodes tend to have most of their connections with weakly connected nodes (see Q for details). 
Disassortativity and other network properties are often used to judge stability of financial networks. 

There has been an explosion in empirical interbank network studies in the last years thanks largely to 
the introduction of electronic settlement systems. One of the first, reported in Q, examines the Austrian 
interbank market, which involves about 900 participating banks. The data are drawn from the Austrian 
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bank balance sheet database (MAUS) and the major loan register (GKE) eontaining all high-value interbank 
loans above €0.36 x 10^; smaller loans are estimated by means of local entropy maximisation. The authors 
construct a network representation of interbank payments for ten quarterly periods from 1999 to 2003. They 
find that the network exhibits small-world properties and is characterised by a power-law distribution of 
degrees. Specifically, the degree distribution is approximated by a power law with the exponent —2.01 for 
degrees > 40. This result, albeit with different exponents, holds for the in- and out-degree distributions too 
(the exponent is —3.1 for out-degrees and —1.7 for in-degrees). A recent study of transactional data from the 
Austrian real-time interbank settlement system (ARTIS) reported in Q demonstrates a strong dependence 
of network topology on the time-scales of observation, with power-law tails exhibiting steeper slopes when 
long time-scales are considered. 

The network structure of transactions between Japanese banks, logged by the Bank of Japan Financial 
Network system (BOJ-NET), is analysed in [l^]. The authors consider several monthly intervals of data 
from June 2001 to December 2002 and construct monthly networks of interbank links corresponding to 21 
transactions or more, i.e. one or more transaction per business day on average. Truncating in this way 
eliminates about 200 out of 546 banks from the network. The resulting monthly networks have a low 
connectivity of 3% and a scale- free cumulative distribution of degrees with the exponent —1.1. 

More than half a million overnight loans from the Italian electronic broker market for interbank deposits 



(e-MID), covering the period from 1999 to 2002, are analysed in [llj. There are about 140 banks in the 
network, connected by about 200 links. The degree distribution is found to exhibit fat tails with power-law 
exponent 2.3 (2.7 for in-degrees and 2.15 for out-degrees), the network is disassortative, with smaller banks 
staying on its periphery. In a related paper 12| , the authors make use of the same dataset to uncover liquidity 
management strategies of the participating banks, given the reserve requirement of 2% on the 23rd of each 
month imposed by the central bank. Signed trading volumes are used as a proxy for the liquidity strategies 
and their correlations are analysed. Two distinct communities supporting the dichotomy in strategy are 
identified by plotting the correlation matrix as a graph. The two communities are mainly composed of large 
and small banks respectively. On average, small banks serve as lenders and large banks as borrowers, but 
the strategies reversed in July 2001, when target interest rates in the Euro area stopped rising and started 
to decrease. The authors also note that some mostly small banks tend to maintain their reserves through 
the maintenance period. The evolution of the network structure over the monthly maintenance period is 
examined in [isj . 

A study of the topology of the Fedwire network, a real-time gross settlement (RTGS) system operated 
by the Federal Reserve System in the USA, is reported in [31 ■ The study covers 62 days in the 1st quarter 
of 2004, during which time Fedwire comprised more than 7500 participants and settled 3.45 x 10^ payments 
daily with total value $1.3 trillion. It reveals that Fedwire is a small- world network with low connectivity 
(0.3%), moderate reciprocity (22%), and a densely connected sub-network of 25 banks responsible for the 
majority of payments. Both in- and out-degree distributions follow a power law for degrees > 10 (exponent 
2.15 for in-degrees and 2.11 for out-degrees). The network is disassortative, with the correlation of out- 
degrees equal to —0.31. The topology of overnight loans in the federal funds market in the USA is examined 



15|, using a large dataset spanning 2415 days from 1999 to 2006. It is revealed that the overnight loans 



form a small- world network, which is sparse (connectivity 0.7%), disassortative (assortativity ranging from 
—0.06 to —0.28), and has low reciprocity of 6%. The reciprocity changes slowly with time and appears to 
follow the target interest rate over the period of several years. A power law is the best fit for the in-degree 
distribution, but the fit is only good for a limited range of degrees. A negative binomial distribution, which 
requires two parameters rather than one for a power law, fits the out-degree distribution best. 

A comprehensive survey of studies of interbank networks is given in [l6| . The number of interbank 
markets being analysed continues to increase. For example, a study of the interbank exposures in Brazil 
for the period from 2004 to 2006 was reported in [13]. A topological analysis of money market flows logged 
in the Danish large- value payment system (Kronos) in 2006 was reported in [T^j, where customer-driven 
transactions are compared with the bank-driven ones. Empirical network studies have been used to guide the 
development of a network model of the interbank market based on the interbank credit lending relationships 

0- . . . . 

Establishing basic topological features of interbank networks is essential for understanding these complex 
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systems. Fundamentally, however, interbank money markets are flow networks, in which links between the 
nodes correspond to monetary flows. The dynamics of such flows has not been examined in depth in previous 
studies, which mostly viewed interbank networks as static or slowly varying. But the underlying flows are 
highly dynamic and complex. Moreover, monetary flows are inhomogeneous; loan flows are fundamentally 
different from the flows of other payments. Payments by the banks' customers and the banks themselves 
cause imbalances in the exchange settlement accounts of the banks. For some banks, the incoming flows 
exceed the outgoing flows on any given day; for other banks, the reverse is true. Banks with excess reserves 
lend them in the overnight money market to banks with depleted reserves. This creates interesting dynamics: 
payment flows cause imbalances, which in turn drive compensating flows of loans. Understanding this 
dynamic relationship is needed for advancing our ability to model interbank markets effectively. 

In this paper, our objective is to define empirically the dynamics of interbank monetary flows. Unlike 
most studies cited above, we aim to uncover the fundamental causal relationship between the flows of 
overnight loans and other payments. We choose to specialise in the Australian interbank market, where 
we have privileged access to a high-quality dataset provided by the Reserve Bank of Australia (RBA). Our 
dataset consists of transactions settled in the period from 19 to 23 February 2007 in the Australian interbank 
market. We separate overnight loans and other payments (which we call nonloans) using a standard matching 
procedure. The loan and nonloan transactions settled on a given day form the flow networks, which are 
the main target of our statistical analysis. We compare the topology and variation of the loan and nonloan 
networks and reveal the causal mechanism that ties them together. We investigate the dynamical stability 
of the system by testing how individual flows vary from day to day. Basic network properties such as the 
degree distribution and assortativity are examined as well. 



2. Data 

High-value transactions between Australian banks are settled via the Reserve Bank Information and 
Transfer System (RITS) operated by the RBA since 1998 on an RTGS basis [l^. The transactions are 
settled continuously throughout the day by crediting and debiting the exchange settlement accounts held 
by the RBA on behalf of the participating banks. The banks' exchange settlement accounts at the RBA 
are continuously monitored to ensure liquidity, with provisions for intra-day borrowing via the intra-day 
liquidity facility provided to the qualifying banks by the RBA. This obviates the need for a monthly reserve 



cycle of the sort maintained by Italian banks as discussed in [13| . The RITS is used as a feeder system 
for transactions orig inating from SWIFT^ and Austraclear for executing foreign exchange and securities 
transactions respectively. The member banks can also enter transactions directly into RITS. The switch to 
real-time settlement in 1998 was an important reform which protects the payment system against systemic 
risk, since transactions can only be settled if the paying banks possess sufficient funds in their exchange 
settlement accounts. At present, about 3.2 x 10^ transactions are settled per day, with total value around 
A$168 billion. 

The data comprise all interbank transfers processed on an RTGS basis by the RBA during the week of 
19 February 2007. During this period, 55 banks participated in the RITS including the RBA. The dataset 
includes transfers between the banks and the RBA, such as RBA's intra-day repurchase agreements and 
money market operations. The real bank names are obfuscated (replaced with labels from A to BP) for 
privacy reasons, but the obfuscated labels are consistent over the week. The transactions are grouped into 
separate days, but the time stamp of each transaction is removed. 

During the week in question, around 2.5 x 10^ transactions were settled per day, with the total value of 
all transactions rising above A$2 x 10^^ on Tuesday and Thursday. The number of transactions (volumcQ) 
and the total value (the combined dollar amount of all transactions) for each day are given in Table [TJ 



^Society for Worldwide Interbank Financial Telecommunication 

■^The term "volume" is sometimes used to refer to the combined dollar amount of transactions. In this paper, we only use 
the term "volume" to refer to the number of transactions and "total value" to refer to the combined dollar amount. This usage 
follows the one adopted by the RBA |20|| . 
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Date 


Volume 


Value (A$10'J) 


19-02-2007 


19425 


82.2506 


20-02-2007 


27164 


206.1023 


21-02-2007 


24436 


161.9733 


22-02-2007 


25721 


212.1350 


23-02-2007 


26332 


184.9202 



Table 1: The number of transactions (volume) and their total value (in units of AS 10^) for each day. 



Date Component 1 Component 2 





(u) 




P 


{u) 




P 


19-02-2007 


4.00 


1.12 


0.81 


6.68 


0.68 


0.19 


20-02-2007 


3.55 


0.72 


0.43 


5.73 


1.49 


0.57 


21-02-2007 


3.66 


0.86 


0.55 


5.86 


1.43 


0.45 


22-02-2007 


3.87 


1.01 


0.68 


6.42 


1.07 


0.32 


23-02-2007 


3.82 


0.87 


0.61 


6.12 


1.19 


0.39 



Table 2: Mean (u), variance a'^, and mixing proportion P of the Gaussian mixture components shown in 
Figure [1] (u = log^o where v is value). 



Figure [T] shows the distribution of transaction values on a logarithmic scale. Local peaks in the distribution 
correspond to round values. The most pronounced peak occurs at A$10^. 

In terms of the number of transactions, the distribution consists of two approximately log-normal com- 
ponents, with lower-value transactions being slightly more numerous. The standard entropy maximisation 
algorithm for a Gaussian mixture model with two components [21j produces a satisfactory fit with the pa- 
rameters indicated in Table [51 The lower- and higher- value components are typically centred around ASlO"* 
and A$10^ respectively. The high-value component is small on Monday (19-02-2007) but increases notice- 
ably on subsequent days, while the low-value component diminishes. By value, however, the distribution is 
clearly dominated by transactions above A$10^, with the highest contribution from around A$2 x 10^. 



3. Overnight loans 

The target interest rate of the RBA during the week of our sample was rt = 6.25% per annum. If 
the target rate is known, it is easy to extract the overnight loans from the data by identifying reversing 
transactions on consecutive days. A hypothetical interest rate can be computed for each reversing transaction 
and compared with the target rate. For instance, suppose a transaction of value vi from bank A to bank B on 
day 1 reverses with value V2, from bank B to bank A, on day 2. These transactions are candidates for the first 
and second legs of an overnight loan from A to B. The hypothetical interest rate for this pair of transactions 
is given by = 100% x 365 x {v2 — vi)/vi] note that the quoted target rate is per annum. Since large banks 
participate in many reversing transactions that can qualify as loans, we consider all possible hypothetical 
pairs and prefer the one that gives rh closest to the target rate. The algorithm for loan extraction is applied 
from Monday to Thursday; loans issued on Friday cannot be processed since the next day is not available. 



A similar procedure was pioneered by Furfine [22| ; see also [23 1 . 

The application of the above algorithm results in the scatter diagram shown in Figure [5] There is a 
clearly visible concentration of the reversing transaction pairs in the region v > 2x 10^ and \rt —rh\< 0.5% 
(red box). We identify these pairs as overnight loans. Contamination from nonloan transaction pairs that 
accidentally give a hypothetical rate close to the target rate is insignificant. By examining the adjacent 
regions of the diagram, i.e. v > 2 x 10^ and rji outside of the red box, we estimate the contamination to be 
less than 2% (corresponding to < 5 erroneous identifications per day) . It is also possible that some genuine 
loans fall outside our selection criteria. However, it is unlikely that overnight interest rates are very different 
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Figure 1: The distribution of transaction values v (in Australian dollars) on a logarithmic scale, with bin 
size Alogj^gW = 0.1; the vertical axis is the number of transactions per bin. Components of the Gaussian 
mixture model are indicated by the dashed curves; the solid curve is the sum of the two components. The 
dotted histogram shows the relative contribution of transactions at different values to the total value (to 
compute the dotted histogram we multiply the number of transactions in a bin by their value). 



from the target rate; and the lower-value transactions (below A$10^), even if they are real loans, contribute 
negligibly to the total value. 

We identify 897 overnight loans over the four days. A daily breakdown is given in Table [31 Here and 
below, we refer to the first leg of the overnight loans as simply loans and to all other transactions as nonloans. 
The loans constitute less than 1% of all transactions by number and up to 9% by value (cf. Tables [T] and [3]) . 
The distribution of loan values and interest rates is shown in Figures l3al and [3bl The interest rate distribution 
peaks at the target rate 6.25%. The mean rate is within one basis point (0.01%) of the target rate, while 
the standard deviation is about 0.07%. The average interest rate increases slightly with increasing value of 
the loan; a least-squares fit yields tk = 6.248 -f 0.010 logio(u/A$10^). 

The same technique can be used to extract two-day and longer-term loans (up to four-day loans for our 
sample of five consecutive days). Using the same selection criteria as for the overnight loans, our algorithm 
detects 27, 67, and 24 two-day loans, with total values A$1.3, A$2.2, and A$1.4 billion, on Monday, Tuesday, 
and Wednesday, respectively. The total value of the two-day loans is 1.5%, 1.0%, and 0.9% of the total 
transaction values on these days respectively. 



4. Nonloans 

We display the distributions of the incoming and outgoing nonloan transactions, for which the bank is 
the destination and the source respectively, for the six largest banks in Figure 3) The distributions are 
similar to the total distribution shown in Figure [U with the notable exception of BA (see below). There is 
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Figure 2: Hypothetical interest rate rh versus value of the first leg of the transaction pairs detected by our 
algorithm, with no restrictions on value or interest rate. The dotted rectangle contains the transactions that 
we identify as overnight loans. The least-squares fit is shown with a solid red line. 



Date 


Volume 


Value (ASlO*^) 


Loan fraction 


19-02-2007 


185 


7.50 


9.12% 


20-02-2007 


221 


9.18 


4.45% 


21-02-2007 


226 


11.08 


6.84% 


22-02-2007 


265 


14.93 


7.04% 



Table 3: Statistics of the overnight loans identified by our algorithm: the number of loans (volume), the 
total value of the first leg of the loans (in units of A$10^), and the fraction of the total value of the loans 
(first legs only) with respect to the total value of all transactions on a given date. 



also an unusually large number of A$106 and A$400 transactions from W to T on Monday. Note that the 
daily imbalance for each bank is mostly determined by the highest value transactions; large discrepancies 
between incoming and outgoing transactions at lower values arc less relevant. 

The distribution for BA is clearly bimodal; it contains an unusually high proportion of transactions 
greater than A$10^. Moreover, below A$10^, incoming transactions typically outnumber outgoing ones by 
a large amount. BA is also involved in many high value transactions that reverse on the same day. These 
transactions probably correspond to the central bank's repurchase agreements, which facilitate intra-day 
liquidity of the banks 2^ . 

The banks shown in Figure|3]are also the largest in term of the number of transactions, with the exception 
of BA. The rank order by the number of transactions matches that by value. For D, which is the largest, 
the number of nonloan transactions reaches 48043 over the week. By the number of transactions, the order 
of the top twelve banks is D, BP, AV, T, W, AH, AF, U, AP, BI, BA, P. By value, the order is D, BP, 
AV, BA, T, W, BG, U, A, AH, AB, BM. The situation is similar when considering the overnight loans. By 
value, AV, D, BP, and T dominate. For these four banks, weekly total loans range from A$11.5 to A$18 
billion and number from 254 to 399. For the other banks the total loan value is less than A$3 billion. 

In view of the discussion above, it is noteworthy that Australia's retail banking system is dominated 
by four big banks (ANZ, CBA, NAB, and WBC jl that in February 2007 accounted for 65% of total resi- 
dent assets, according to statistics published by Australian Prudential Regulation Authority (APRA); see 



http://www.apra.gov.au for details. The resident assets of the big four exceeded A$225 billion each, well 
above the next largest retail bank, St George Bank Limitec(f| (A$93 billion). The distinction between the big 



^Australia and New Zealand Banking Group, Commonwealth Bank of Australia, National Australia Bank, and Westpac 
Banking Corporation. 

*In December 2008, St George Bank became a subsidiary of Westpac Banking Corporation. 
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Figure 3: (a) The distribution of loan values w on a logarithmic scale. The vertical axis is the number of 
loans per bin for bin size Alogj^g v = 0.25. The dotted line is the same distribution multiplied by the value 
corresponding to each bin (in arbitrary units). The date of the first leg of the loans is indicated. |(b)| The 
distribution of loan interest rates r^. The vertical axis is the number of loans per bin for bin size Ar^ — 0.01. 
The date of the first leg of the loans is indicated. The mean and standard deviation are 6.25% and 0.08% 
on Monday (19-02-2007), and 6.26% and 0.07% on the other days. 



four and the rest of the banks in terms of cash and liquid assets at the time was less clear, with Macquarie 
Bank Limited in third position with A$8 billion. According to APRA, cash and liquid assets of the big four 
and Macquarie Bank Limited accounted for 56% of the total. 



5. Loan and nonloan imbalances 

In order to maintain liquidity in their exchange settlement accounts, banks ensure that incoming and 
outgoing transactions roughly balance. However, they do not control most routine transfers, which are 
initiated by account holders. Therefore, the imbalances arise. On any given day, the nonloan imbalance of 
bank i is given by 

j fc i k 

where {wfc(i,.7)}fc is a list of values of individual nonloan transaction from bank i to bank j, settled on the 
day. The nonloan imbalances are subsequently compensated by overnight loans traded on the interbank 
money market. The loan imbalances are defined in the same way using transactions corresponding to the 
first leg of the overnight loans. Note that we do not distinguish between the loans initiated by the banks 
themselves and those initiated by various institutional and corporate customers. For instance, if the funds 
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Figure 4: The distribution of uonloan transaction values of the six largest banks for Monday through 
Thursday (from left to right); the banks are selected by the combined value of incoming and outgoing 
transactions over the entire week. Black and red histograms correspond to incoming (bank is the destination) 
and outgoing (bank is the source) transactions; red histograms are filled in to improve visibility. The banks' 
anonymous labels, the combined daily value of the incoming and outgoing transactions, and the daily 
imbalance (incoming minus outgoing) are quoted at the top left of each panel (in units of A$10^). The 
horizontal axis is the logarithm of value in A$. 



of a corporate customer are depleted, this customer may borrow overnight to replenish the funds. In this 
case, the overnight loan is initiated by an account holder, who generally has no knowledge of the bank's net 
position. Nevertheless, the actions of this account holder in acquiring a loan reduce the bank's imbalance, 
provided that the customer deposits the loan in an account with the same bank. 

The loan and nonloan imbalances for the six largest banks are given in Tabled The data generally comply 
with our assumption that the overnight loans compensate the daily imbalances of the nonloan transactions. 
The most obvious exception is for BA on Thursday (22-02-2007), where a large negative nonloan imbalance 
is accompanied by a sizable loan imbalance that is also negative. Taking all the banks together, there is a 
strong anti-correlation between loan and nonloan imbalances on most days. We see this clearly in Figure [SJ 
The Pearson correlation coefficients for Monday through Thursday are —0.93, —0.88, —0.95, —0.36. It is 
striking to observe that many points fall close to the perfect anti-correlation line. The anti-correlation is 
weaker on Thursday (crosses in Figure [5]), mostly due to BA and AV. 

A correlation also exists between the absolute values of loan imbalances and the nonloan total values 
(incoming plus outgoing nonloan transactions); the Pearson coefficients are 0.74, 0.75, 0.66, 0.77 for Monday 
through Thursday. This confirms the intuitive expectation that larger banks tolerate larger loan imbalances. 
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Av log^o V 

Figure 5: Left: loan imbalance AZ vs nonloan imbalance Av for individual banks and days of the week (in 
units of A$10^). Right: the absolute value of loan imbalance \Al\ vs nonloan total value (incoming plus 
outgoing transactions) for individual banks and days of the week. Thursday data are marked with crosses. 





19-02-2007 


20-02-2007 


21-02-2007 


22-02-2007 




nonloans 


loans 


nonloans 


loans 


nonloans 


loans 


nonloans 


loans 


D 


-0.51 


+0.12 


-0.28 


+0.12 


-0.76 


+0.45 


-1.25 


+ 1.44 


BP 


+2.08 


-1.64 


+0.80 


-0.59 


+1.38 


-0.85 


+0.16 


+0.82 


AV 


-0.32 


-0.17 


+1.39 


-0.79 


+0.55 


-0.65 


+1.08 


+0.25 


BA 


+0.03 


-0.19 


-0.10 


-0.31 


-0.32 


-0.05 


-1.53 


-0.64 


T 


-0.76 


+ 1.10 


-0.75 


+0.68 


-0.62 


+0.64 


-0.21 


+0.27 


W 


-0.09 


+0.07 


+0.08 


+0.26 


-0.36 


+0.41 


+0.36 


-0.37 



Table 4: Loan and nonloan imbalances for the six largest banks (in units of A$10^). 



6. Flow variability 

For each individual source and destination, we define the nonloan flow as the totality of all nonloan 
transactions from the given source to the given destination on any given day. The value of the flow is the 
sum of the nonloan transaction values and the direction is from the source to the destination. On any given 
day, the value of the flow from bank i to bank j is defined by 

Vfiow{i,j) = ^Vk{i,j), (2) 

k 

where {vk{i, j)}k is a list of values of individual nonloan transaction from i to j on the day. For example, 
all nonloan transactions from D to AV on Monday form a nonloan flow from D to AV on that day. The 
nonloan transactions in the opposite direction, from AV to D, form another flow. A flow has zero value if 
the number of transactions is zero. Typically, for any two large banks there are two nonloan flows between 
them. The loan flows are computed in a similar fashion. 

6. 1 . Nonloan flows 

There are 55 banks in the network, resulting in A^aow = 2970 possible flows. The actual number of flows 
is much smaller. The typical number of nonloan flows is ^ 800 on each day (the actual numbers are 804, 
791, 784, 797). Even though the number of nonloan flows docs not change signiflcantly from day to day, 
we find that only about 80% of these fiows persist for two days or more. The other 20% are replaced by 
different flows, i.e. with a different source and/or destination, on the following day. Structurally speaking, 
the network of nonloan flows changes by 20% from day to day. However, persistent flows carry more than 
96% of the total value. 

Even when the flow is present on both days, its value is rarely the same. Given that 80% of the network 
is structurally stable from day to day, we assess variability of the network by considering persistent flows 
and their values on consecutive days. Figure [5] shows the pairs of persistent flow values for Monday and 



logiQ V on 19-02-2007 log^o v on 20-02-2007 log^p f on 21-02-2007 

Figure 6: Nonloan flow value pairs on one day (horizontal axis) and the next (vertical axis). Only flows 
present on both days arc considered. Flows that do not change lie on the diagonal (red dotted line). The 
solid line is the weighted orthogonal least squares fit to the scatter diagram; the weights have been defined 
to emphasize points corresponding to large flows. 



Tuesday, Tuesday and Wednesday, and Wednesday and Thursday. If the flow values were the same, the 
points in Figure [S] would lie on the diagonals. We observe that the values of some flows vary signiflcantly, 
especially when comparing Monday and Tuesday. Moreover, there is a notable systematic increase in value 
of the flows from Monday to Tuesday by a factor of several, which is not observed on the other days. For 
each pair of days shown in Figure |6j we compute the Pearson correlation coefficient, which gives 0.53 for 
Monday and Tuesday, 0.70 for Tuesday and Wednesday, and 0.68 for Wednesday and Thursday. 

To characterize the difference between the flows on different days more precisely, we compute the Eu- 
clidean distance between normalised flows on consecutive days. We reorder the adjacency matrix {wflow(*, j)}ij 
of the flow network on day d as an A'flow-dimensional vector v^; representing a list of all flows on day d 
{d — 1,2,..., 5). For each pair of consecutive days we compute the Euclidean distance between normalized 
vectors Vd/lvd] and Vd+i/|vd+i|, which gives 0.62, 0.50, 0.50 for all flows and 0.61, 0.49, 0.49 for persis- 
tent flows (the latter are computed by setting non-persistent flows to zero on both days). Since the flow 
vectors are normalized, these quantities measure random flow discrepancies while systematic deviation such 
as between the flows on Monday and Tuesday are ignored. For two vectors of random values uniformly 
distributed in interval (0, 1), the expected Euclidean distance is 0.71 and the standard deviation is 0.02 for 
the estimated number of persistent nonloan flows of 640. So the observed variability of the nonloan flows is 
smaller than what one might expect if the flow values were random. 

6. 2. Loan flows 

Variability of the loan flows is equally strong. The number of loan flows varies from 69 to 83 (actual 
numbers are 69, 75, 77, 83). Only about 50% of these flows are common for any two consecutive days. 
Moreover, persistent flows carry only about 65% of the total value of the loan flows on any given day, cf. 
80% of nonloan flows. For persistent loan flows, the Pearson correlation coefficients are 0.63, 0.90, and 0.76 
for the consecutive pairs of days starting with Monday and Tuesday. The correlation is generally similar to 
that of the nonloan flows, with the notable exception of the loan flows on Tuesday and Wednesday, when 
the sub-network of persistent loan flows appears to be more stable. 

The Euclidean distances between the normalized loan flows for each pair of consecutive days are 0.85, 
0.68, 0.73 for all flows and 0.63, 0.44, and 0.44 for persistent flows. For two vectors of random values 
uniformly distributed in interval (0, 1), the expected Euclidean distance is 0.7 and the standard deviation 
is 0.1 for the estimated number of persistent loan flows of 40. So the observed variability of the persistent 
loan flows is much smaller than what one might expect if the flow values were random. 
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Figure 7: As for Figure [6] but for loan flows. 

6.3. Relation between nonloan and loan flows 

Some loan flows do not have corresponding nonloan flows between the same nodes on the same day. These 
flows carry about 14% of loan value on Monday, and about 7% on Tuesday through Thursday. Nonloan 
flows that have corresponding loan flows account for 35% to 48% of all nonloan flows by value, even though 
the number of these flows is less than 10% of the total. 

To improve the statistics, we aggregate the flows on all four days. Figure [8] shows nonloan and cor- 
responding loan flow values. We fail to find any significant correlation between loan and nonloan flows 
(Pearson coefficient is 0.3). The correlation improves if we restrict the loan flows to those consisting of 
three transactions or more; such flows mostly correspond to large persistent flows. In this case the Pearson 
coefficient increases to 0.6; banks that sustain large nonloan flows can also sustain large loan flows, even 
though the loan flows on average are an order of magnitude lower than the corresponding nonloan flows. 
The lack of correlation when all loans are aggregated is due to the presence of many large loans that are 
not accompanied by large nonloan transactions, and vice versa. 

7. Net flows 

The net flow between any two banks is defined as the difference of the opposing flows between these banks. 
The value of the net flow equals the absolute value of the difference between the values of the opposing flows. 
The direction of the net flow is determined by the sign of the difference. If vnowihj) > Vfiov,{j,i), the net 
flow value from i to j is given by 

Wnct(i, j) = Wflow(i, - fflow(j,«)- (3) 

For instance, if the flow from D to AV is larger than the flow in the opposite direction, then the net flow is 
from D to AV. 

7.1. General properties 

The distributions of net loan and nonloan flow values are presented in Figure The parameters of the 
associated Gaussian mixture models are quoted in Table [5l The distribution of net nonloan flow values 
has the same general features as the distribution of the individual transactions. However, unlike individual 
transactions, net flow values below ASlO'' are rare; net flows around A$10^ are more prominent. 

There are on average around 470 net nonloan flows each day. Among these, roughly 110 consist of a 
single transaction and 50 consist of two transactions, mostly between small banks. At the other extreme, 
net flows between the largest four banks (D, BP, AV, T) typically have more than 10"^ transactions per day 
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Figure 8: Loan flow values versus nonloan flow values combined over four days. Triangles correspond to 
loan flows with three or more transactions per flow. The solid line is the orthogonal least squares fit to the 
scatter diagram; the weighting is the same as in Figure |6l 



Date Component 1 Component 2 





{u) 




P 


{u) 




P 


19-02-2007 


5.14 


1.88 


0.60 


7.51 


0.36 


0.40 


20-02-2007 


5.70 


2.17 


0.51 


7.82 


0.50 


0.49 


21-02-2007 


5.73 


1.97 


0.52 


7.72 


0.44 


0.48 


22-02-2007 


5.78 


2.06 


0.57 


7.86 


0.45 


0.43 



Table 5: Mean (it), variance ct^, and mixing proportion P of the Gaussian mixture components appearing 
in Figure ini (u = logio v). 

each. Overall, the distribution of the number of transactions per net flow is approximated well by a power 
law with exponent a — —1.0 ± 0.2: 

7V„et(n) cxn", (4) 

where N^^tin) is the number of net nonloan flows that consist of n transactions (n ranges from 1 to more 
than 1000). This is consistent with the findings for Fcdwirc reported in [l5| (see right panel of Fig. 14 in 

There are roughly 60 net loan flows each day. As many as 40 consist of only one transaction. On the 
other hand, a single net loan flow between two large banks may comprise more than 30 individual loans. The 
distribution of the number of transactions per net loan flow is difficult to infer due to poor statistics, but it 
is consistent with a power law with a steeper exponent, —1.4 ± 0.2, than that of the nonloan distribution. 
There are no net loan flows below A$10^ or above A$10^. Comparing net loan and nonloan flows, it is 
obvious that net loan flows cannot compensate each and every net nonloan flow. Not only are there fewer 
net loan flows than nonloan flows, but the total value of the former is much less than the total value of the 
latter. 
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Figure 9: The distribution of values of net nonloan flows (black histogram) on a logarithmic scale with bin 
size Alogj^g V = 0.1. The components of the Gaussian mixture model are indicated with the dashed curves; 
the solid curve is the sum of the two components. Net loan flows are overplotted in red. The vertical axis 
counts the number of net flows per bin. 



Net loan and net nonloan flows arc not correlated; the correlation coefficient is 0.3. Restricting net loan 
flows to those that have three transactions or more does not improve the correlation. If a net loan flow 
between two banks was triggered to a significant degree by the magnitude and the direction of net nonloan 
fiow between these bank, one expects a correlation between net loan and nonloan fiows. Our examination 
shows that in this respect loan flows are decoupled from nonloan flows. The connection between them is 
indirect. Namely, nonloan flows cause an imbalance in the account of each bank, which is subsequently 
compensated by loan flows, which are largely unrelated to the nonloan flows that caused the imbalance. 



7.2. Degree distribution and assortativity 

We define the in-degree of node i as the number of net flows that terminate at i, i.e. the number of net 
flows with destination i, and the out-degree as the number of net flows that originate from i, i.e. the number 
of net flows with source i. The degree distribution of the nonloan networks is shown in Figure [TUal Node 
BA has the highest in-degree of 37 on Monday, but on the other days it drops to 15 on average, while the 
out-degree is 11.75 on average for this node. The highest in-degrees are usually found among the four largest 
banks (D, BP, AV, T); the only exception is Monday, when AF's in-degree of 22 is greater than AV's 21, 
and BA has the highest in-degree. The highest out-degrees are usually achieved by D, BP, AV, T, W, and 
AH; the exceptions are Monday, when D's out-degree of 17 is less than AR's and AP's 18, and Thursday, 
when AV's out-degree of 16 is less than P's 18. 

It is difficult to infer the shape of the degree distribution for individual days due to poor statistics. The 
two-sample Kolmogorov-Smirnov (KS) test does not distinguish between the distributions on different days 
at the 5% significance level. With this in mind, we combine the in- and out-degree data for all four days 
and graph the resulting distributions in Figure [TObl We find that a power law distribution does not provides 
a good fit for either in- or out-degrees. Visually, the distribution is closer to an exponential. However, the 
exponential distribution is rejected by the Anderson-Darling test. 

The degree distribution conceals the fact that fiows originating or terminating in nodes of various degrees 
have different values and therefore provide different contributions to the total value of the net flows. Nodes 
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Figure 10: (a) Degree distribution of the net nonloan flow networks (for convenience, in-degrees are positive 



and out-degrees are negative) . The total value of the net flows corresponding to the specific degrees is shown 
with red dots (the log of value in A$10^ is indicated on the right vertical axis). |(b)| Degree distribution of the 
net nonloan flows when the degree data for all four days are aggregated (in-degrees are circles; out-degrees 
are triangles). 



with lower degrees are numerous, but the flows they sustain are typically smaller than those carried by a few 
high-degree nodes. In particular, for the nonloan flows, nodes with in-degree d < 10 are numerous, ranging 
from 35 to 37, but their outgoing net flows carry about 20% of the value on average. On the other hand, 
nodes with d > 17 are rare, but their flows carry 50% of the value. The same effect is observed for the 
out-degrees. 

The degree distribution of the network of net loan flows is shown in Figure [TTal fwe ignore the nodes that 
have zero in- and out- degrees over four days). Similarly to nonloan flows, the KS test does not distinguish 
between the distributions on different days at the 5% signiflcance level. The combined distribution is shown 
in Figure lllbl 

To probe assortativity of the net flow networks, we compute the in-assortativity deflned in (25| as the 
Pearson correlation coefficient between the in-degrees of sources and destinations of the net flows (out- 
assortativity is computed similarly using the out-degrees). The net nonloan flow network is disassortative, 
with in-assortativity of —0.39, —0.37, —0.38, —0.37 and out-assortativity of —0.35, -0.38, —0.39, —0.37 on 
Monday, Tuesday, Wednesday, and Thursday, respectively. The net loan flow network is less disassortative; 
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Figure 11: (a) Same as Figure ITUal but for the net foan flow networks. |(b)|Same as Figure ITObl but for the 



net loan flows. 



the in-assortativity is -0.16, -0.26, -0.18, -0.19 and the out-assortativity is -0.03, -0.10, 0.02, -0.20 for 
the same sequence of days. In biological networks, the tendency of out-assortativity to be more assortative 
than in-assortativity has been noted in 25 1. 



7.3. Topology of the net flows 

Given the source and destination of each net flow, we can construct a network representation of the 
net flows. An example of the network of net nonloan flows is shown in Figure [121 The size of the nodes 
and the thickness of the edges are proportional to the net imbalances and net flow values respectively (on 
a logarithmic scale). We use the Fruchterman-Reingold algorithm to position the nodes [26|; the most 
connected nodes arc placed in the centre, and the least connected nodes are moved to the periphery. The 
core of the network is dominated by the four banks with the largest total value and the largest number 
of transactions: D, BP, AV, and T. The other big banks, such as AF, AH, and W, also sit near the core. 
It is interesting to note the presence of several poorly connected nodes (Q, V, BF, and especially X) that 
participate in large incoming and outgoing flows, which produce only negligible imbalances in the banks 
themselves. 

The sub-network consisting of D, BP, AV, BA, T, W, U, A, AH, AF, AP, and P is fully connected on 
all five days, i.e. every node is connected to every other node. The sub-network of D, AV, and BP is fully 
connected, even if we restrict the net flows to values above A$10^. 
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Figure 12: Network of net nonloan flows on Tuesday, 20-02-2007. White (grey) nodes represent negative 
(positive) imbalances. The bank labels are indicated for each node. The size of the nodes and the thickness 
of the edges are proportional to the logarithm of value of the imbalances and the net flows respectively. 

In Figure [T^l the flows between the largest nodes are difficult to discern visually, because the nodes are 
placed too close to each other in the image. We therefore employ the following procedure to simplify the 
network. We consider the fully connected sub-network of twelve nodes, plus node BG, and combine all other 
nodes into a new node called "others" in such a way that the net flows are preserved (BG is included because 
it usually participates in large flows and is connected to almost every node in the complete sub-network). 
The result of this procedure applied to the daily nonloan networks is presented in Figures [13aH13d1 For these 
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plots, we employ the weighted Fruchterman-Reingold algorithm, which positions the nodes with large flows 
between them close to each other. The imbalances shown in Figure I13bl are the same as those of the full 
network in Figure [T21 The daily networks of net loan flows for the same nodes are shown in Figures ri4aH 14dl 
We observe that the largest flows on Monday (19-02-2007) were significantly lower than the flows on 
the subsequent days. The largest nodes (D, BP, AV, T, W) are always placed close to the center of the 
network, because they participate in the largest flows. The topology of the flows is complex and difficult 
to disentangle, even if one concentrates on the largest flows (above A$5 x 10®). For instance, on Monday, 
probably the simplest day, the flow of nonloans is generally from BG to "others" to D to BP. There are also 
sizable flows from T to AV and from AV to "others" and BP. However, lower value flows (below A$5 x 10®) 
cannot be neglected completely because they are numerous and may contribute significantly to the imbalance 
of a given node. 

Nodes D, T, BP, AV, and W form a complete sub-network of net loan flows on Monday, Tuesday, and 
Wednesday. This sub-network is almost complete on Thursday too, except for the missing link between BP 
and W. The appearance of the net loan network is different from that of the nonloan network, since the 
same nodes participate in only a few loan flows. Therefore, the position of a node in the network image 
is strongly influenced by the number of connections of that node. Some of the poorly connected nodes are 
placed at the periphery despite the fact that they possess large flows. The four largest nodes (D, T, BP, 
AV) are always positioned at the center of the network. 

l.J^. Network variability 

The net nonloan flow network is extremely volatile in terms of flow value and direction. For example, a 
A$10^ flow from D to BP on Monday transforms into a A$3.2 x 10^ flow in the same direction on Tuesday, 
only to be replaced by a A$6.3 x 10® flow in the opposite direction on Wednesday, which diminishes further 
to A$2.5 X 10^ on Thursday. Nodes T and BP display a similar pattern of reversing flows between Tuesday 
and Wednesday. On the other hand, the net flow between T and AV maintains the same direction, but the 
flow value is strongly fluctuating. In particular, a moderate A$4.8 x 10® flow on Monday rises to A$1.9 x 10^ 
on Tuesday, then falls sharply to A$2 x 10® on Wednesday and again rises to A$2.2 x 10^ on Thursday. 

Considering any three nodes, we observe that circular and transitive flows are present on most days, the 
latter being more common. The most obvious example is a circular flow between D, T, and BP on Thursday 
and a transitive flow involving BG, T, and AV on the same day. The circular flows are unstable in the sense 
that they do not persist over two days or more. 

The net loan flow network exhibits similar characteristics. Few net loan flows persist over the four days. 
For example, the flow from AV to T has the same direction and is similar in value on all four days. Circular 
loan flows are also present, as the flow between AV, T, and BP on Thursday demonstrates. 

8. Conclusions 

In this paper, we study the properties of the transactional flows between Australian banks participating 
in RITS. The value distribution of transactions is approximated well by a mixture of two log- normal compo- 
nents, possibly reflecting the different nature of transactions originating from SWIFT and Austraclear. For 
the largest banks, the value distributions of incoming and outgoing transactions are similar. On the other 
hand, the central bank displays a high asymmetry between the incoming and outgoing transactions, with 
the former clearly dominating the latter for transactions below A$10^. 

Using a matching algorithm for reversing transactions, we successfully separate transactions into loans 
and nonloans. For overnight loans, we estimate the identification rate at 98%. The mean derived interest 
rate is within 0.01% of the central banks' target rate of 6.25%, while the standard deviation is about 0.07%. 
We find a strong anti-correlation between loan and nonloan imbalances (Pearson coefficient is about 0.9 on 
most days). A likely explanation is that nonloan flows create surpluses in some banks. The banks lend the 
surplus to banks in deficit, creating loan flows that counteract the imbalances due to the nonloan flows. 
Hence, loan and nonloan imbalances of individual banks are roughly equal in value and opposite in sign on 
any given day. 
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(c) 21-02-2007 



(d) 22-02-2007 



Figure 13: Networks of daily net nonloan flows for D, AV, BP, T, W, BA, AH, AF, U, AP, P, A, BG. All 
the other nodes and the flows to and from them are combined in a single new node called "others". The 
size of the nodes and the thickness of the edges are proportional to the logarithm of value of the imbalances 
and the net flows respectively. The value of the flows and the imbalances can be gauged by referencing a 
network shown in the middle, where the values of the flows are indicated in units of A$l billion. 
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(c) 21-02-2007 (d) 22-02-2007 



Figure 14: Networks of daily net loan flows. The same nodes as in Figures ri3aH13dl are used. The scale of 
the loan flows, the imbalances, and the positions of the nodes are the same as those used for the nonloan 
flows in Figures fT3aHl3dl to simplify visual comparison. 
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The flow networks arc structurally variable, with 20% of nonloan flows and 50% of loan flows replaced 
every day. Values of persistent flows, which maintain the same source and destination over at least two 
consecutive days, vary significantly from day to day. Some flow values change by several orders of magnitude. 
Persistent flows increase in value several-fold between Monday and Tuesday. Individual flow values can 
change by several orders of magnitude on the following day. Overall, there is a reasonable correlation 
between the flow values on consecutive days (Pearson coefficient is 0.65 for nonloans and 0.76 for loans 
on average). We also find that larger banks tend to sustain larger loan flows, in accord with the intuitive 
expectations. However, there is no correlation between loan and nonloan flows. 

We examine visually the topology of the net loan and nonloan flow networks. The centre of both networks 
is dominated by the big four banks. Twelve banks form a complete nonloan sub-network, in which each 
bank is connected to every other bank in the sub-network. The three largest banks form a complete sub- 
network even if the net flows are restricted to values above ASlO'^. Our examination reveals that the network 
topology of net flows is complicated, with even the largest flows varying greatly in value and direction on 
different days. 

Our findings suggest a number of avenues for future research on interbank networks. Firstly, the relation- 
ships we uncovered can be used to constrain analytical models and numerical simulations of interbank flows 
in financial networks. In particular, our explanation of the link between the loan and nonloan imbalances 
needs to be tested in numerical simulations. Secondly, it is necessary to analyse interbank markets in other 
countries to establish what elements of our results are signatures of general dynamics and what aspects are 
specific to the epoch and location of this study. Even when high quality data are available, most previous 
studies concentrate on analysing static topological properties of the networks or their slow change over time. 
The internal dynamics of monetary flows in interbank networks has been largely ignored. Importantly, 
one must ask whether the strong anti-correlation between loan and nonloan imbalances is characteristic of 
RTGS systems whose institutional setup resembles the Australian one or whether it is a general feature. 
For instance, in Italy a reserve requirement of 2% must be observed on the 23rd of each month, which may 
encourage strong deviations between loan and nonloan imbalances on the other days. 

Acknowledgement 

We thank the Reserve Bank of Australia for supplying the data. AS acknowledges generous flnancial 
support from the Portland House Foundation. 

References 

[1] E. Kolaczyk, Statistical analysis of network data: methods and models, Springer Verlag, 2009. 
[2] M. Jackson, Social and economic networks, Princeton Univ Pr, 2008. 

[3] G. Caldarelli, Scale-Free Networks: Complex webs in nature and technology, Oxford University Press, USA, 2007. 

[4] F. Schweitzer, G. Fagiolo, D. Sornette, F. Vega-Redondo, D. White, Economic networks: What do we know and what do 

wc need to know?, Advances in Complex Systems 12 (4-5) (2009) 407-422. 
[5] A. Haldane, R. May, Systemic risk in banking ecosystems, Nature 469 (7330) (2011) 351-355. 
[6] R. May, S. Levin, G. Sugihara, Ecology for bankers, Nature 451 (21) (2008) 893-895. 
[7] M. Newman, Mixing patterns in networks, Physical Review E 67 (2) (2003) 26126. 

[8] M. Boss, H. Elsinger, M. Summer, S. Thurner, Network topology of the interbank market. Quantitative Finance 4 (6) 
(2004) 677-684. 

[9] F. Kyriakopoulos, S. Thurner, C. Puhr, S. Schmitz, Network and eigenvalue analysis of financial transaction networks, 
The European Physical Journal B-Condcnsed Matter and Complex Systems 71 (4) (2009) 523-531. 

[10] H. Inaoka, T. Ninomiya, K. Taniguchi, T. Shimizu, H. Takayasu, Fractal network derived from banking transaction-an 
analysis of network structures formed by financial institutions. Bank of Japan Working Papers 4. 

[11] G. Dc Masi, G. fori, G. Caldarelli, Fitness model for the Italian interbank money market, Physical Review E 74 (6) (2006) 
066112. 

[12] G. lori, R. Reno, G. De Masi, G. Caldarelli, Trading strategies in the Italian interbank market, Physica A: Statistical 

Mechanics and its Applications 376 (2007) 467-479. 
[13] G. lori, G. De Masi, O. Precup, G. Gabbi, G. Caldarelli, A network analysis of the Italian overnight money market, 

Journal of Economic Dynamics and Control 32 (1) (2008) 259-278. 
[14] K. Soramaki, M. L. Bech, J. Arnold, R. J. Glass, W. Beyeler, The topology of interbank payment flows, Physica A: 

Statistical Mechanics and its Applications 379 (1) (2007) 317-333. 



20 



[15] M. L. Bcch, E. Atalay, The topology of the federal funds market, Physica A: Statistical Mechanics and its Applications 
389 (22) (2010) 5223-5246. 

[16] K. Imakubo, Y. Socjima, The transaction network in Japanese interbank money markets, Monetary and Economic Studies 
28 (2010) 107-150. 

[17] D. Cajueiro, B. Tabak, The role of banks in the Brazilian interbank market: Does bank type matter?, Physica A; Statistical 

Mechanics and its Applications 387 (27) (2008) 6825-6836. 
[18] K. R0rdam, M. Bech, The topology of Danish interbank money flows, FRU Working Papers. 

[19] S. Li, J. He, Y. Zhuang, A network model of the interbank market, Physica A: Statistical Mechanics and its Applications. 
[20] P. Gallagher, J. Gauntlett, D. Sunner, Real-time Gross Settlement in Australia, Reserve Bank of Australia Bulletin 

September (2010) 61. 
[21] G. McLachlan, D. Peel, Finite mixture models, Wiley-Intersciencc, 2000. 

[22] C. Furflne, Interbank Exposures: Quantifying the Risk of Contagion., Journal of Money, Credit & Banking 35 (1) (2003) 
111-129. 

[23] A. Ashcraft, D. Duffie, Systemic illiquidity in the federal funds market. The American economic review 97 (2) (2007) 
221-225. 

[24] F. Campbell, Reserve Bank Domestic Operations under RTGS, Reserve Bank of Australia Bulletin November (1998) 54. 
[25] M. Piraveenan, M. Prokopenko, A. Zomaya, Assortative Mixing in Directed Biological Networks, IEEE IEEE/ACM 

Transactions on Computational Biology and Bioinformatics. 
[26] T. Fruchterman, E. Reingold, Graph drawing by force-directed placement, Software - Practice and Experience 21 (11) 

(1991) 1129-1164. 



21 



